Measurements of the j*p ~ ¥ A Reaction At Low Q 2 : Probing the Mesonic Contribution 



S. Stave/'LjN. Sparveris, 2 <[i] M. O. Distlcr, 3 I. Nakagawa, 1 ' 4 '!! P. Achcnbach, 3 C. Ayerbe 
Gayoso, 3 D. Baumann, 3 J. Bcrnaucr, 3 A. M. Bernstein, x 'l R. Bohm, 3 D. Bosnar, 5 T. Botto, 1 
A. Christopoulou, 2 D. Dale, 4 M. Ding, 3 L. Doria, 3 J. Friedrich, 3 A. Karabarbounis, 2 M. Makek, 5 

H. Merkel, 3 U. Miiller, 3 R. Neuhausen, 3 L. Nungesser, 3 C.N. Papanicolas, 2 A. Piegsa, 3 
J. Pochodzalla, 3 M. Potokar, 6 M. Seimetz, 3 '0 S. Sirca, 6 S. Stiliaris, 2 Th. Walcher, 3 and M. Weis 3 

(Al Collaboration) 

department of Physics, Laboratory for Nuclear Science and Bates Linear Accelerator 
Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA 
2 Institute of Accelerating Systems and Applications and Department of Physics, University of Athens, Athens, Greece 
3 Institute fur Kernphysik, Johannes Gutenberg- Universitat Mainz, D-55099 Mainz, Germany 
^Department of Physics and Astronomy, University of Kentucky, Lexington, Kentucky 40206 USA 
5 Department of Physics, University of Zagreb, Croatia 
6 Institute Jozef Stefan, University of Ljubljana, Ljubljana, Slovenia 
(Dated: July 2, 2008) 

The determination of non-spherical angular momentum amplitudes in nucleons at long ranges 
(low Q 2 ), was accomplished through the p(e,e'p)ir° reaction in the A region at Q 2 = 0.060, 0.127, 
and 0.200 (GeV/c) 2 at the Mainz Microtron (MAMI) with an accuracy of 3%. The results for 
the dominant transition magnetic dipole amplitude and the quadrupole to dipole ratios have been 
obtained with an estimated model uncertainty which is approximately the same as the experimental 
uncertainty. Lattice and effective field theory predictions agree with our data within the relatively 
large estimated theoretical uncertainties. Phenomenological models are in good agreement with 
experiment when the resonant amplitudes are adjusted to the data. To check reaction model cal- 
culations additional data were taken for center of mass energies below resonance and for the o TL i 
structure function. These results confirm the dominance, and general Q 2 variation, of the pionic 
contribution at large distances. 
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I. INTRODUCTION 

Experimental confirmation of the presence of non- 
spherical hadron amplitudes (i.e. d states in quark mod- 
els or p wave 7T-N states) is fundamental and has been the 
subject of intense experimental and theoretical interest 
(for reviews see (T), H, H, H, 5f ) . This effort has focused on 
the measurement of the electric and Coulomb quadrupole 
amplitudes (E2, C2) in the predominantly Ml (magnetic 
dipole-quark spin flip) ^y*N — > A transition. 

The present low Q 2 experiments add important data to 
determine the physical basis of long range nucleon and 
A non-spherical amplitudes. This is the region where 
pionic effects are predicted to be dominant and appre- 
ciably changing. The experiment was carried out at the 
Mainz Microtron to measure cross sections and extract 
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the resonant multipoles at Q 2 = 0.060, 0.127, and 0.200 
(GeV/c) 2 . The Q 2 = 0.060 (GcV/c) 2 point is the low- 
est Q 2 value probed to date in modern electroproduc- 
tion experiments. The Q 2 = 0.200 (GeV/c) 2 point tests 
the Q 2 variation and provides a valuable overlap with 
newly obtained Jefferson Lab data 0. The Q 2 = 0.127 
(GeV/c) 2 point tested the background amplitudes and 
acted as a comparison of results from Mainz and Bates. 
Aspects of this work arc given in @, [H . This work in- 
cludes more of the details and is a complete account of 
those data and includes previously unpublished data as 
well. 

The present measurements fill in an important gap in 
the coverage of the Q 2 evolution between the photon 
point (Q 2 = 0) 0, Mi 

and previously published electro- 
production experiments at JLab [ll|, 0~1, G3 f° r Q 2 from 
0.4 to 6.0 (GcV/c) 2 , with the exception of good coverage 
at Q 2 = 0.127 (GeV/c) 2 at Bates 0, El O E3 and 
Mainz [jj, [H, 0] at Q 2 = 0.127, 0.200 (GeV/c) 2 that 
have been published. 

Since the proton has spin 1/2, no quadrupole moment 
can be measured. However, the A has spin 3/2 so the 
ry* N — > A reaction can be studied for quadrupole am- 
plitudes in the nucleon and A. Due to spin and par- 
ity conservation in the 7*A r (J 7r = 1/2+) -> A(J ff = 
3/2 + ) reaction, only three multipoles can contribute to 
the transition: the magnetic dipole (Ml), the electric 
quadrupole (E2), and the Coulomb quadrupole (C2) 
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photon absorption multipoles. The corresponding res- 
onant pion production multipoles arc M l+ 

3 /2 

S 1+ . The relative quadrupole to dipole 

EMR=Re(£^{ 2 /Mi 3 { 2 ) and CMR^Ref^/Mff ). In 
the quark model, the non-spherical amplitudes in the nu- 
cleon and A are caused by the non-central, tensor interac- 
tion between quarks [2l], [22j| ■ However, the magnitudes 
of this effect for the predicted E2 and C2 amplitudes [23| 
are at least an order of magnitude too small to explain 
the experimental results (see Fig. [TCI below) and even the 
dominant Ml matrix clement is ~ 30% low 0, [23J. A 
likely cause of these dynamical shortcomings is that the 
quark model does not respect chiral symmetry, whose 
spontaneous breaking leads to strong emission of virtual 
pions (Nambu-Goldstone Bosons) [2]. These couple to 
nucleons as a ■ p where a is the nucleon spin, and p is 
the pion momentum. The coupling is strong in the p 
wave and mixes in non-zero angular momentum compo- 
nents. Based on this, it is physically reasonable to expect 
that the pionic contributions increase the Ml and domi- 
nate the E2 and C2 transition matrix elements in the low 
Q 2 (large distance) domain. This was first indicated by 
adding pionic effects to quark models 0, 125I I2II . sub- 
sequently in pion cloud model calculations [271 . |28| | , and 
recently demonstrated in effective field theory (chiral) 
calculations 123. [3011. 



II. EQUIPMENT 

The p(e, e'p)ir° measurements were performed using 
the Al spectrometers at the Mainz Microtron [3l|. Elec- 
trons were detected in Spectrometer A which used two 
pairs of vertical drift chambers for track reconstruction 
and two layers of scintillator detectors for timing infor- 
mation and particle identification. The protons were de- 
tected in Spectrometer B which has a detector package 
similar to Spectrometer A. Spectrometer B also has the 
ability to measure at up to 10° out-of-plane in the lab. 
Due to the Lorentz boost, this corresponds to a signifi- 
cantly larger value in the center of mass frame. The mo- 
mentum resolution of the spectrometers is 0.01% and the 
angular resolution at the target is 3 mrad [3l|. Details 
about the spectrometers are available in 3l| . The MAMI 
B accelerator delivered a longitudinally polarized, contin- 
uous, electron beam up to 855 MeV. Beam polarization 
was measured periodically with a M0ller polarimeter [32j 
to be ~ 75%. The beam with average current of up to 
25 fiA was scattered from a liquid hydrogen cryogenic 
target. The beam energy has an absolute uncertainty of 
±160 keV and a spread of 30 keV (FWHM) [3l|. The 
effects of these uncertainties and the various kinematic 
cuts were studied to estimate an overall systematic un- 
certainty (see Table HI)) for the cross sections of 3 to 4%. 
This was tested with elastic electron-proton scattering 
and the data agree with a fit to the world data [33[ at 
the 3% level. In addition, a third spectrometer (Spec- 



trometer C) was used throughout the experiment as a 
luminosity monitor. 



III. EXPERIMENTAL METHODOLOGY 

The five-fold differential cross section for the 
p(e, e'p)TT° reaction is written as five two-fold differential 
cross sections with an explicit 0* dependence as [34| 
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where 0* 9 is the pion center of mass azimuthal angle 
with respect to the electron scattering plane, h is the he- 
licity of the electron beam, p e is the polarization of the 
electron beam, vlt = ^/2e(l + e), vlt> = \/2e(l — e), e 
is the transverse polarization of the virtual photon, and T 
is the virtual photon flux. The virtual photon differential 
cross sections (<tt, &l,&lt, &tt, glt 1 ) ar c all functions of 
the center of mass energy W, the four momentum trans- 
fer squared Q 2 , and the pion center of mass polar angle 
8* q (measured from the momentum transfer direction). 
They are bilinear combinations of the multipoles [34| . 

The extraction of the cross sections was performed us- 
ing three sequential measurements. For the helicity in- 
dependent cross sections there are three cross sections 
to extract: ctq = <jt + (e was not varied so the 
two cross sections cannot be separated), cttt and ctlt- 
The three two-fold differential cross sections can be ex- 
tracted algebraically by measuring the five-fold differen- 
tial cross section at the same center-of-mass energy W, 
four-momentum transfer squared Q 2 , and proton center- 
of-mass polar angle 9* but different values of the proton 
azimuthal angle 4>* q . (The proton and pion are back-to- 
back in the center of mass frame leading to the follow- 
ing relations between the angles: 0* — 180° — 6** 9 and 
<f> Pq = 180° + <t% q .) The sequential kinematic settings 
then keep W and Q 2 constant by keeping the electron 
arm (Spectrometer A) unchanged and the proton arm 
was moved so that 91 q remained the same and 4>* q was 
changed. Using Eq.[T| the three two- fold differential cross 
sections (ctq, <jtt, &lt) can then be found algebraically 
from the three measured five-fold cross sections and the 
(j)* q angles at which they were measured. The fifth two- 
fold differential cross section ctlt' was measured by re- 
versing the helicity of the longitudinally polarized elec- 
tron beam at non-zero (out-of-plane) 4>* q angles, ctlt' is 
sensitive to the background terms and provides another 
test of the reaction calculations. 

Figure [1] shows the kinematic overlap for the sequen- 
tial (j>* pq settings at Q 2 = 0.060 (GeV/c) 2 . The W over- 
lap is approximately 40 MeV, the AQ 2 « 0.04 GeV 2 /c 2 , 
A0* q w 10°, and A<j)* pq w 40°. For Q 2 = 0.200 (GcV/c) 2 , 
the overlap region is slightly larger due to the larger 
Lorentz boost but the shapes are qualitatively similar. 
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FIG. 1: (Color online) Plot of the overlap of the sequential settings for Q 2 = 0.060 (GeV/c) 2 ,W = 1221 MeV, 6* pq = 24°. 
Medium gray corresponds to <f>p q = 0°, light gray to (j>p q = 90° and dark gray to 4>vq = 180°. The amount of overlap for 
Q 2 = 0.20 (GeV/c) 2 is qualitatively similar. 



Studies of the extraction process showed that the 
smallest uncertainties and most sensitivity were achieved 
when the three <j)* q measurements were as far apart as 
possible. However, at the larger 9* q angles, not all cj>* q 
values can be reached because of the 10° out-of-plane an- 
gle constraint of Spectrometer B. Therefore, for each 9* q 
setting, the maximally out-of-plane settings were used. 
However, each kinematic setting was carefully chosen in 
order to minimize the uncertainties in the algebraic cross 
section extraction process. 



The kinematics for all of the setups are shown in detail 
in Table Q] The data presented in this work were taken 
during two run periods in 2003. The first period was in 
April and measured the mostly non-parallel cross sections 
for Q 2 = 0.060 and 0.200 (GeV/c) 2 in addition to an 
extension of the Q 2 = 0.127 (GeV/c) 2 data set. The 
October period was used to measure the W scans at Q 2 = 
0.060 and 0.200 (GeV/c) 2 and the low W background 
terms. In addition to the Q 2 = 0.060,0.127, and 0.200 
(GeV/c) 2 measurements, some cross check measurements 
were made at Q 2 = 0.127 (GeV/c) 2 which overlapped 
with existing data from Bates. 



TABLE I: Kinematic values for W, Q 2 , proton center-of-mass 
polar angle 9* q , proton azimuthal angle (f>* q and the initial 
electron beam energy Ef,e am . See Tables IIVI and IVHII for 
detailed settings for the W scans. 
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IV. DATA ANALYSIS 
A. Phase Space Acceptance and Simulation 

The phase space acceptance in the spectrometers is 
in a multi-dimensional space and has a complex shape 
(see Fig. [TJ. One challenge was defining the phase space 
acceptance in a similar manner across all the kinematic 
settings with the phase space varying by a large amount 
across the spectrometer acceptance. One solution is to 
have a very small acceptance which will limit the varia- 
tions but also limit the statistics. Too large of a phase 
space acceptance leads to large systematic errors as the 
variation in phase space is too large for the simulation to 
reliably calculate. However, a compromise can be found 
by settling in a region where the combination of the sta- 
tistical and systematic errors is a minimum. This is il- 
lustrated in Fig. [2] where the normalized cross section is 
plotted against the size of the kinematic cut. More de- 
tails of the calculation are presented below but the effect 
of too small of a phase space acceptance (large statistical 
errors on the left side of Fig. \2§ can be seen as well as the 
effect of too large of a phase space acceptance (large sys- 
tematic errors seen as deviations from the central value 
on the right side of Fig. [5]). 

To ensure uniformity of the phase space selection 
across the varying kinematics, a unique solution was 
found which non-arbitrarily defined the edges of the ac- 
ceptance. The maximum allowed phase space region was 
found by locating the half-maximum points in the dis- 
tributions of the variables upon which the cross section 
depends: W, Q 2 , 9* q , 4>* q . Symmetry around the central 
kinematics was also enforced so that neither side of the 
phase space was weighted too heavily. Once the max- 
imum acceptance regions were defined, then fractional 
widths of those regions were used to study the behavior 
of the extracted cross section. Those studies, detailed be- 
low, were then used to define the final phase space region 
used for extracting the cross sections. 

SimulH — h [HI is the software which was employed to 
calculate the multi-dimensional phase space. Simul+-|- 
also simulates the collimators inside the spectrometers 
as part of calculating the phase space. After the subtrac- 
tion of background events (see next subsection) , the spec- 
trometer acceptance was limited in software to the cen- 
tral region of the spectrometers to keep edge effects out of 
the analysis. The details are in (36|. In addition to pre- 
cise spectrometer properties and collimators, Simul+-|- 
also calculates energy loss and the radiative corrections 
in the same way as for the data. Each simulated event 
contains the proper weighting for radiative corrections, 
the virtual photon flux T and the lab to center-of-mass 
Jacobian. The simulated events undergo kinematic selec- 
tion processes identical to those used on the data and can 
then be used to determine the phase space and, finally, a 
cross section. 

Figure [3] shows a comparison of the data for an 
in-plane, forward setup with the results of Simul++ 
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FIG. 2: Variation in the cross section due to changes in cut 
size for Q 2 = 0.060 (GeV/c) 2 , 0* pq = 24°, <f)* vq = 0°. The 
abscissa shows the fractional phase space selection width for 
all the variables mentioned in the text and the ordinate shows 
the cross section normalized to the cross section result with 
the 0.50 fractional phase space selection width. The bars are 
centered on the cut corrected result. The dark bar is the sta- 
tistical uncertainty and the light bar shows the total uncer- 
tainty with the systematic uncertainty added in quadrature. 



weighted by the MAID 2003 phenomcnological model [37[ 
cross section and plotted against the four physics vari- 
ables upon which the cross section depends: W, Q 2 , 9* q , 
4>* q . As is clear in the figure, there is very good agree- 
ment for all the variables across the acceptance. A fifth 
variable, z, was also examined closely because it affects 
the size of the spectrometer acceptance, z is the vertex 
position determined by Spectrometer B which has better 
vertex resolution than Spectrometer A. The real edges of 
the z distribution are not as sharp as in the simulation, 
but extensive studies showed that avoiding those regions 
in z yielded reliable cross section results. Other setups 
have similarly good agreement between data and simu- 
lation. In addition to good agreement in the previously 
listed variables, there is also acceptable agreement on the 
shape and location of the missing mass peak as shown in 
Fig. [U The differences between simulation and data for 
the missing mass do not cause appreciable uncertainties 
and the level of agreement is sufficient for this analysis. 

To investigate the effects of different sized phase space 
acceptance regions on the extracted cross section, sev- 
eral types of studies were performed. In all of them, 
the maximum phase space selection width was defined 
by the half-maximum points as mentioned earlier. Then, 
only the fractional width of the phase space region rela- 
tive to the maximum width was varied. In this manner, 
the phase space selection was consistent across the many 
kinematic settings. 

Fine scans were made for each kinematic setting by si- 
multaneously varying the fractional width of the phase 
space selection of all the physics variables (W , Q 2 , 6* q , 
(j)* q ) and the vertex position, z. (The cross section does 
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FIG. 3: (Color online) Comparison of the relative shapes in phase space for the data and the simulation weighted by MAID 
2003 [33 for the primary variables used in the cut. The results are for an in-plane, forward angle setup at Q 2 = 0.060 
(GeV/c) 2 and show a good level of agreement. 
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FIG. 4: (Color online) Comparison of the missing mass from 
the data (black crosses) and the simulation (solid line) for 
the Q 2 = 0.060 (GeV/c) 2 , W = 1221 MeV, 0* pq = 36°, 
4>p q — 180° setting. 



not depend on the vertex position but the shape of the 
vertex distribution did change from setup to setup re- 
quiring a similar definition of the cut.) Figure [2] shows 
the scan for a non-parallel, forward angle, Q 2 = 0.060 
(GeV/c) 2 setting. The abscissa shows the fractional 
phase space selection width for all the variables men- 
tioned and the ordinate shows the cross section normal- 
ized to the cross section result with the 0.50 fractional 



phase space selection width. The variation of the cross 
section ratio with changing cut fractions seen in Fig. [2] 
was representative of the variation seen in the other kine- 
matic settings. As mentioned before, Figure [2] shows that 
the extracted cross section gets more stable with smaller 
selection regions but the statistical uncertainties neces- 
sarily get larger. Small statistical uncertainties are pos- 
sible with larger cut fractions but then the systematic 
errors suffer. What is not shown here is that the helic- 
ity dependent cross sections have the most stable results 
for a fractional phase space selection width of 0.75. The 
fraction of 0.75 was then chosen as a compromise in or- 
der to have stable results across all kinematics with small 
statistical uncertainties. 

In most settings, like that shown in Fig. [3J the cross 
section for a fractional phase space selection width of 
0.75 is slightly lower than the 0.50 result. To correct for 
this, a phase space correction factor was determined by 
averaging the results of these selection scans over all the 
similar kinematic settings. These corrections are all on 
the order of 2 to 3% with a 1% systematic uncertainty. 
(Some of the backward angle settings had flat phase space 
selection width scan results and did not require a phase 
space correction factor.) For a comparison of the rela- 
tive uncertainties, see the light and dark bars in Fig. [5] 
The light bar shows the statistical uncertainty for the 
final cross section result with a fractional phase space se- 
lection width of 0.75 and the dark bar shows the total 
uncertainty including the appropriate systematic uncer- 
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taintics. (Since this is a comparison of only one kinematic 
setting, it is not appropriate to include any uncertainties 
from quantities which vary statistically or systematically 
from setup to setup like the luminosity.) The cross sec- 
tion ratios with fractional phase space selection widths 
of 0.50 and 0.75 agree with each other within statistical 
and systematic uncertainty and are also both stable. 

This phase space selection procedure was used for the 
analysis of the Q 2 = 0.060 (GeV/c) 2 data. Similar cuts 
were used in the Q 2 — 0.200 (GeV/c) 2 analysis and their 
stability was verified. Any small differences in the cuts 
lead to only small differences in the cross sections and 
are not significant. 
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B. Elimination of Background Counts 

During the pion production runs, there were two types 
of backgrounds: 7r~ / fj,~ background and general acciden- 
tal background. The ir~ / [i~ background was removed by 
making a two-dimensional selection in missing mass ver- 
sus coincidence timing space. The ir~ background 
region was very clearly separated from ir° events of in- 
terest. Their identity was confirmed with a Chercnkov 
counter in Spectrometer A which was only present during 
the first running period. However, the two-dimensional 
selection in missing mass versus coincidence timing space 
was found to be just as effective at removing the j 
background. 

After the 7r~//i~ cut, an accidental subtraction was 
applied using accidentals from both sides of the coin- 
cidence peak to determine the background counts per 
channel. Figure [5] shows the coincidence peak on top of 
the accidental background with the light gray region indi- 
cating the average background level seen in the two side 
regions. The accidental subtraction removes about 6% 
to 20% of the events in the coincidence peak depending 
on the kinematics and is the largest of the background 
subtractions. 

After both the pion and accidental background sub- 
tractions, the data consist of only coincidence events and 
a cross section can be extracted which is not contami- 
nated with background events. 



C. Luminosity 

The luminosity is calculated based upon the total cur- 
rent measured by the Forster probe, a pair of toroidal 
coils which surround the beam and measure the current 
induced by the beam [38|]. The Forster probe is located 
in the third stage of the microtron which can recircu- 
late the beam up to 90 times. Therefore, the current of 
the recirculated beam in the third stage can be up to 90 
times larger than the beam on target. A measurement 
there leads to a much more precise determination of the 
beam current. 



FIG. 5: Coincidence timing of the proton and electron. The 
dark gray areas indicate the background subtraction region 
and the light gray under the peak indicates the size of the 
background. The full-width at half-maximum of the peak is 
0.9 ns. 



The luminosity can then be calculated given the beam 
current, target length, and target density (from pressure 
and temperature) . To prevent local boiling of the hydro- 
gen target, the electron beam is rasterized or wobbled 
across the target in a rectangular pattern. During the 
April run, the beam was also placed off-center to ensure 
a path to the out-of-plane Spectrometer B that was free 
of obstructions. The flat plate above the target extended 
out and would have been in the path of the out-going 
protons if the beam were not shifted down and to the 
right. This offset in the beam position decreased the ef- 
fective target length by less than 1.5% and the effect was 
taken into account by the simulation. 

The normal operating pressure for the target is 2.1 bar. 
With a normal temperature of 22 K, this leads to an un- 
dercooling of 1 K. This temperature buffer allows for a 
certain amount of local heating without the target start- 
ing to boil. However, both pion production run periods 
experienced lower target pressure which led to less un- 
dercooling. Instead of 1 K undercooling, the experiment 
operated closer to 0.6 K undercooling. 

The singles rates in Spectrometer A were used to study 
the effect of the beam current on the luminosity. (Spec- 
trometer B was rarely in the same place from one setup 
to another but Spectrometer A was returned to the same 
location repeatedly.) By plotting the singles rate in A 
versus beam current, any target boiling effect should be 
visible. Figure [5] shows the A singles rate divided by 
the beam current for all the Q 2 — 0.060 (GeV/c) 2 runs. 
Also plotted are the average and the RMS deviation of 
all of the data. Notice that almost all of the points are 
consistent with a horizontal line which indicates no beam 
current dependent luminosity change. Therefore, the low 
Q 2 runs were below the boiling threshold and do not need 
any correction. 

However, other data were taken with higher beam cur- 
rents, specifically the parallel pion production cross sec- 
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FIG. 6: (Color online) The singles rate in Spectrometer 
A divided by the beam current for all the Q 2 = 0.060 
(GeV/c) 2 runs in the April beam time period. The uncer- 
tainties shown are statistical only. The lines show the average 
and RMS deviation of all the data points. Most of the data 
are within the uncertainties. Even over a range of 3 /J.A, there 
is no large target boiling effect. 



tion comparison with Bates data and the Q 2 — 0.200 
(GeV/c) 2 data. It is possible that these runs were taken 
above the boiling threshold. To test this, the singles rates 
from Spectrometer A divided by the current was plotted 
and a line was fit to the data. During the experiment, the 
effect of the beam current on the luminosity was explic- 
itly checked for one setting where data were taken at 25 
fiA for 2.5 hours and 12.5 fiA for 5 hours. The results of 
the fit to the singles data and the beam current study in- 
dicate a current dependent effect for beam currents above 
12.5 fiA. A luminosity correction factor and uncertainty 
of (3 ± 1)% were adopted which are consistent with all 
the available data. More details are presented in [36[ . 

The conclusion from the luminosity studies is that the 
Q 2 = 0.06 (GeV/c) 2 data are unaffected by beam target 
heating and a (3 ± 1)% correction is sufficient to account 
for the effect in the remaining data. 



D. Extraction at Central Kinematic Values 

The analysis procedure yields a cross section which has 
been averaged over the multi-dimensional phase space 
while the theoretical models provide predicted values at 
points in that phase space. To compare the averaged 
cross section to theory, a kinematic translation procedure 
is applied to the data. This is also known as bin centering 
corrections [H| or transport. The goal of the procedure 
is to find the correction factor which will convert the 
cross section which has been averaged over phase space 
to the cross section evaluated at the central kinematic 
values of the phase space. The kinematic correction fac- 
tor is found by averaging the model predictions over the 
same volume in phase space as the data. That value is 



then divided by the model prediction at the center of the 
phase space. The inverse of that ratio is the correction 
factor. This technique does not rely upon the absolute 
size of the theory but merely requires that the theory 
have the same shape throughout the same phase space 
as the data. Corrections are typically 2 to 3% indicating 
that the cross section tends to vary smoothly and fairly 
symmetrically through the phase space. A small (0.5%) 
systematic uncertainty is introduced with this method 
which was estimated by performing the translation with 
several models and taking the RMS deviation of the re- 
sults. 

This method of translation was tested by varying the 
size of the phase space selection region and checking for 
convergence to the point cross section. Smaller cuts lead 
to larger statistical uncertainties but the tests showed 
that the results were stable and converged within the 
uncertainties. 



E. Absolute Cross Section Verification 

In order to determine stability over time and the 
proper normalization, the elastic reaction p(e, e'p) was 
measured throughout the experiment. As during the pion 
production runs, Spectrometer A was used to detect elec- 
trons and Spectrometer B for protons. The measurement 
uncertainties are dominated by the systematics estimated 
at approximately 4%. The results are stable over time 
and are consistent, within systematic uncertainties, with 
the 1996 dispersion-theoretical analysis fit to the world 
clastic data |33|. The 2004 dispersion analysis 1401 1 and 




other fits to the elastic scattering data [41 
were examined and there is only a small amount of spread 
between the fits and they agree at (98.5 ± 1.5)% of the 
1996 dispersion fit. A more recent dispersion analysis [45[ 
is slightly lower (about 95% of the 1996 fit) but agrees 
with the other fits and the data within the systematic 
errors. 

The conclusion from the coincidence elastic analysis is 
that the measured cross sections are stable over time and 
agree well with previous elastic results. This indicates 
stability in the luminosity, target density, and beam po- 
sition. It also indicates that the spectrometers can be 
placed reliably (typically 0.6 mm and 0.1 mrad [1H) and 
that the central momenta are well known. 

For another check of the absolute cross sections, the 
parallel pion production cross section at W — 1232 MeV 
and Q 2 = 0.127 (GeV/c) 2 was taken during both run 
times. This measurement had been carried out previ- 
ously at Bates [H, [13, IH, 53, 51] . The beam energy for 
each of these past experiments was slightly different and 
so the e factor is slightly different. This can be corrected 
for using the ratio of &t to ctl from a model. Using 
MAID 2003, the correction factor is about 1% and is 
even smaller for other models. Figure [7] shows all of the 
parallel cross section comparisons for the previous Bates 
data and the current experiment. There is a reasonable 
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FIG. 7: (Color online) Comparison of the parallel cross section 
from previous experiments at Bates (o: Mertz-Vellidis (MV) 
e = 0.614 and Sparveris (Spar.) e = 0.768 0) and from 
the current experiment (•) all measured at or converted to 
e = 0.707. The uncertainties are the statistical and systematic 
uncertainties added in quadrature. The / lines are the average 
of the central values and the uncertainties for Bates. The \ 
lines are the values for Mainz. The overlap region is easily 
seen. 



overlap region since the systematic uncertainties are ac- 
counted for in the plot. The Mainz results arc stable over 
time from April 2003 to October 2003. Another item to 
consider is that the variation in the Bates measurements 
is about 4% and the difference from Mainz to the lowest 
Bates point is about the same. The conclusion drawn is 
that the current measurement agrees with previous mea- 
surements within the systematic uncertainties. 



F. Systematic Uncertainties 

As mentioned above, the uncertainties from the kine- 
matic translation procedure can be estimated by using 
various models and looking at the RMS deviation and, 
for most settings, the effect is less than 0.5%. The one 
exception is a 1.3% effect due to worse phase space over- 
lap in the Q 2 = 0.06 (GeV/c) 2 , 6 pq , = 29.6° setting. 
This was caused by unforeseen difficulties in placing the 
spectrometers. 

Table HTl summarizes the remaining uncertainties. The 
luminosity uncertainty comes from a 1% uncertainty in 
the target length and a 1% uncertainty in the density. 
Those estimates have been conservatively added linearly. 
However, the stability of the elastic cross sections indi- 
cates that this systematic uncertainty should affect all 
runs in the same way. The detector inefficiency cor- 
rection was estimated in previous works and is quoted 
here [H, |49[ . The dead time correction factor was cal- 
culated using vetoed and unvetoed scalers and is based 
upon counting statistics. 

The phase space cut uncertainties were found by vary- 



Uncertainty 


Size [%] 


Luminosity 


2 


Detector inefficiency correction 


1 


Dead time correction uncertainty per setup 


< 0.5 


Phase space cut uncertainty 


1.5 - 2.4 


Model uncertainty 


0.4,1.3 


Beam current luminosity correction 


0,1 


Momentum and angular resolution 


1 


Beam position 


1 


Total in quadrature 


3.3 - 3.7 


Beam polarization 


1.2 



TABLE II: Summary of systematic and model uncertainties. 



ing the size of the kinematic phase space cuts. The large, 
in-plane angle settings had very little difference, but for 
the rest of the settings, the difference was typically 2 
to 3%. The systematic uncertainty in phase space cut 
uncertainties was estimated to be the average of the un- 
certainties in the ratios of the small and large cut regions. 
The systematic uncertainty in the cut correction is be- 
tween 1.5% and 2.4%. 

The model uncertainty in kinematic translation has 
already been detailed as has the beam current luminos- 
ity correction uncertainty. Note that the beam current 
related luminosity correction is not applied for beam cur- 
rents less than 12.5 fiA and when it is applied, has a 1% 
uncertainty. 

To see the effect of the spectrometer angular and mo- 
mentum resolution, the central momentum and angle set- 
tings for the spectrometers were shifted in the simulation 
and the shifted simulation results were used to extract 
cross sections. Using various combinations of the resolu- 
tions and for several, representative setups, the resolution 
uncertainty was estimated at 1%. The spectrometer po- 
sitioning uncertainties of 0.6 mm and 0.1 mrad [3l| are 
much smaller than the resolution uncertainties and so do 
not affect the results. The beam position can also af- 
fect the cross section. A study showed that this effect is 
about 1%. 

To summarize, there are several corrections applied to 
the data (luminosity, phase space, kinematic translation) 
but they have all been studied in detail and their contri- 
butions are all well determined. The total systematic 
uncertainties are in the range of 3 to 4% and agree very 
well with the estimates based upon comparisons with the 
world clastic cross sections. 



V. EXPERIMENTAL RESULTS 

The methods described in the previous sections were 
applied to the data at Q 2 = 0.060,0.127 and 0.200 
(GeV/c) 2 and the results are given here and in the Ap- 
pendix in tabular form. There are two types of cross 
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sections presented. One is the five-fold differential cross 
section which is dependent on W, Q 2 , 9* q , and (j>* q and is 
measured directly by the spectrometers. The other type 
is the two-fold differential cross section which is 4>* q in- 
dependent and must be extracted from the five-fold cross 
sections using Eq.[TJ Both types of cross sections are used 
to aid in comparison with theory and for fitting purposes. 



A. Near resonance: Q 2 = 0.060,0.200 (GeV/c) 2 

The extracted partial cross sections <7o, <Jtt, &lt, and 
o LT > versus 6* q for W = 1221 MeV, Q 2 = 0.060 and W = 
1232 MeV, Q 2 = 0.200 (GeV/c) 2 are plotted in Figs. 
and [9] respectively. These data are compared with the 
chiral effective field theory calculations (EFT) which have 
a few low energy parameters and then rely on theory to 
arrive at results. These calculations have relatively large 
estimated uncertainty bands due to the neglect of higher 
order terms. Within these uncertainties the agreement 
with experiment is good. While these calculations and 
their uncertainty estimates are a great contribution to 
the field, conclusions cannot be drawn unless there are 
further improvements. The precision of the data is such 
that more precise theory is required. The inclusion of 
even higher order terms appears to be necessary. 

The top sections of Figures [8] and [9] also show the pre- 
dictions of four model calculations. The Sato-Lee (SL) 
[27[ and Dubna-Mainz- Taipei (DMT) [H[ models contain 
explicit pion cloud contributions while the MAID [37} 
and SAID [50| calculations are primarily phenomenolog- 
ical. These models have been adjusted by their authors 
to agree with our previous data 0, EE EE E3- For 
Q 2 = 0.060 (GeV/c) 2 , all models agree with the data for 
<jtt- For (To only MAID is not in agreement with the 
data. However, for er^y the dispersion between the mod- 
els and data is greater showing that they have not been 
adjusted to agree with Si+. For Q 2 = 0.200 (GeV/c) 2 , 
the agreement between models and experiment for Si+ 
is even less satisfactory. One item that this indicates is 
that the Q 2 dependence of is not correct in DMT 
and MAID since both models agree well with Si+ data 
at Q 2 = 0.127 (GeV/c) 2 shown in Ref. For a LT >, 

only the SL model agrees with the data at both Q 2 val- 
ues. All of these disagreements show the importance of 
performing measurements at low Q 2 . 

The extraction of the three resonant 7* + p — > A 
amplitudes Mi+, E±+ and Si+ was accomplished by 
adjusting these amplitudes in the four phenomcnolog- 
ical models described above. Following the practice 
of Refs. @, H, EE E3 the model dependent extraction 
from successful phenomenological reaction models al- 
lows for a reliable extraction of the resonant amplitudes. 
The model uncertainty is estimated by the spread of 
the derived values using the various model amplitudes 

[i mm. 

The fitting procedure used in this analysis is described 
in detail in [36[ • Briefly, the procedure takes all the 



background multipoles up to L — 5 from a model and 
varies the amplitude of the resonant, isospin 3/2 multi- 
poles (Aff^ 2 , E^J 2 and Sf+ 2 ) to attain a best fit to data 
at one value of W and Q 2 . By performing the fit in this 
manner, there is not the usual truncation of the fit past 
p waves. However, there is a model dependence since the 
various models differ in the sizes of background terms. 

The fitting of the data started with the helicity inde- 
pendent results, the three 8* q angles with the <fi* q depen- 
dence. Those seven five-fold differential cross section re- 
sults were fit using the three resonant parameter fit with 
the four models. All the fits had x 2 P er degree of freedom 
near one indicating good fits. Correlations between the 
fitting parameters were taken into account in the uncer- 
tainties estimated by the fitting routine [IE [|53| . Figures 
[5] and [5] show the data and the different fitted models. 
Despite different background terms, the four model fits 
converged. It is impressive that the four model curves 
almost fall on top of each other when the three resonant 
7*p — > A amplitudes (Mj 3 ^ 2 , E^ 2 , S\+) are varied to fit 
the data as shown in the lower panel of Figs. [5] and [5] 
In addition, the lower panels show the "spherical" calcu- 
lated curves when the resonant quadrupole amplitudes 

3/2 3/2 

(E 1+ in Co and S 1+ in olt) are set equal to zero. The 
difference between the spherical and full curves shows the 
sensitivity of these cross sections to the quadrupole am- 
plitudes and demonstrates the basis of the present mea- 
surement. The small spread in the spherical curves in- 
dicates their sensitivity to the model dependence of the 
background amplitudes. 

Figure [TU] shows the model convergence for the EMR 
and CMR at Q 2 = 0.060 (GeV/c) 2 in another way. The 
convergence in M\ + was not as significant but the values 
for Mi_|_ have been modified by the various model authors 
in order to fit previous data. 

It is interesting that good fits were achieved in the reso- 
nant region despite the differing model backgrounds. The 
Ml 4- term is dominant but the background multipoles are 
of a similar size to the resonant multipoles. The reason 
that the fitting routine is able to be rather insensitive to 
the backgrounds is due, mostly, to their having a differ- 
ent phase. Near resonance, the 7 = 3/2 resonant mul- 
tipoles are mostly imaginary due to the Fermi- Watson 
theorem [34j |. The Eq + and £0+ are mostly imaginary 
while others are primarily real. Since the M\ + amplitude 
near resonance is almost pure imaginary, the interference 
with mostly real amplitudes is very small. In addition, 
the Eq+ multipole does not differ very much from model 
to model so while it has a large effect, it does not affect 
the resonant fits. The fitting procedure is also insensitive 
to the background amplitudes partly because of their an- 
gular dependence. The primary contributors to the cross 
section near resonance are the resonant Mi+, E\+ and 
S\+ and the background Eq + and Sq+. The multipole 
contributions to the cross section have different angular 
shapes which the fitter can use to separate the compo- 
nents. 
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FIG. 8: (Color online) The measured <ro = ctt + £ijl,<Jtt,o'lt, and a hT t differential cross sections as a function of 0* q at 
W = 1221 MeV and Q 2 = 0.060 (GeV/c) 2 . The ▼ symbols are our data points and include the statistical and systematic 
uncertainties added in quadrature. The top figures (panels a-d) show the data with the EFT predictions [3(| which are 
plotted with their estimated uncertainties. The other curves represent predictions from the MAID 2003 [37I. SL(Sato-Lee) [27j [. 
DMT [2p |. and SAID [5(| models. The bottom figures (panels e-h) show our data with model curves for which the three 

3/23/23/2 3/2^3/2 

resonant multipoles M 1+ , E 1+ , S 1+ are fit to the data. The lines with dots are the fitted models with the E 1+ and S 1+ 
quadrupole terms set to zero and are only plotted for the sensitive observables, ctq and ctlt- 



As mentioned before, in Figures [H] and [SJ the ult 1 re- 
sults are only close for the Sato-Lee model but then those 
cross sections were not included in this fit. The ctlt' cross 
section is sensitive primarily to the background ampli- 
tudes and a resonant fit is not expected to improve the 
agreement. In fact, the fit results were the same, within 
the uncertainties, whether or not the <jlt' data were in- 
cluded. 

Table Mil shows model and chiral EFT predictions 
along with fitted results for the models and the aver- 
ages of those models at both Q 2 values. The table also 



contains three different types of uncertainties: statistical 
(used when fitting the data), systematic, and model. The 
systematic uncertainties are calculated by scaling all of 
the cross sections to the minimum and maximum allowed 
by the uncertainties and refitting. The range of the refit 
values then gives the systematic uncertainty. The sys- 
tematic uncertainty for the EMR and CMR mostly can- 
celled because the quantities are ratios of multipoles and 
so are supressed in Table IIIII However, since M\ ( 2 is 
not a ratio, the systematic uncertainties remained. Fol- 
lowing our previous work [J, [l?], HH, HH , the model un- 
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FIG. 9: (Color online) The measured <ro = ctt + £ijl,<Jtt,o'lt, and a LT i differential cross sections as a function of 0* q at 
W = 1221 MeV and Q 2 = 0.200 (GeV/c) 2 . The ▼ symbols are our data points and include the statistical and systematic 
uncertainties added in quadrature. The top figures (panels a-d) show the data with the EFT predictions [3(J which are 
plotted with their estimated uncertainties. The other curves represent predictions from the MAID 2003 [3?| . SL(Sato-Lee) [27f ] . 
DMT [2q |. and SAID models. The bottom figures (panels e-h) show our data with model curves for which the three 

3/23/23/2 3/2*-(3/2 

resonant multipoles M 1+ , E 1+ , S 1+ are fit to the data. The lines with dots are the fitted models with the E 1+ and S 1+ 
quadrupole terms set to zero and are only plotted for the sensitive observables, <tq and olt- 



certainties were found by taking the root mean square 
deviation of the results using the four models. We be- 
lieve that this is reasonable since the chosen models rep- 
resent state-of-the-art calculations and also a variety of 
different approaches. The final statistical and systematic 
uncertainties are the average over the four models. The 
model uncertainties and experimental uncertainties are 
very similar in size, especially for the EMR and CMR, 
as also seen at Q 2 = 0.127 (GcV/c) 2 in Ref. [13 • There- 
fore, one can conclude that the current experimental limit 
has been reached and further gains can only be achieved 



after improving the models. The effect of background 
amplitudes on the resonant amplitudes was studied and 
determined to have an effect approximately the same size 
as the model to model RMS deviation. This study is de- 
tailed in Refs. [H] and 
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FIG. 10: (Color online) Example of the convergence of the EMR and CMR values with fitting for W = 1232 MeV, Q 2 = 0.06 
(GeV/c) 2 (panels a-b) and W = 1232 MeV, Q 2 = 0.20 (GeV/c) 2 (panels c-d). The uncertainty on the fits is statistical only 
since the systematic uncertainty is very small as mentioned in the text. The left side of each plot shows the original model 
calculations and the right side shows the results after fitting. T he g ray band is the total statistical, systematic and model 
uncertainty added in quadrature. The models are MAID 2003 [13, EJ, DMTpl, [HI, Sato-Lee [13] and SAID Jp. The chiral 
effective field theory predictions of Pascalutsa and Vanderhaeghen (PV) [3^, Ugand Gail and Hemmert (GH) |29(| are included. 
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Q 2 (GeV/c) 2 Model 


EMR (%) 


OMR (%) 


Mil 2 (10- 3 /"V+) 

TTif OrifT 
r 11 Ul lg,. 


0.06 SAID 
SL 
DMT 
MAID 


-2.18 ±0.31 -1.80 
-2.26 ±0.30 -2.98 
-2.11 ±0.28 -2.84 
-2.56 ±0.27 -2.16 


-4.87 ±0.29 -5.30 
-4.46 ±0.25 -3.48 
-4.85 ±0.26 -5.74 
-5.07 ±0.26 -6.51 


40.81 ±0.29 ±0.57 40.72 
40.20 ±0.27 ±0.56 41.28 
40.78 ±0.27 ±0.57 40.81 
39.51 ±0.26 ±0.57 40.53 


Avg. 


-2. 28 ±0.29 ±0.01 ±0.20 


-4.81 ±0.27 ±0.03 ±0.26 


40.33 ±0.27 ±0.57 ±0.61 


GH 
PV 


-2.66 

-2.88 ± 0.70 


-6.06 

-5.85 ± 1.40 


41.15 

39.75 ± 3.87 


Q 2 (GeV/c) 2 Model 


EMR (%) 
Fit Orig. 


CMR (%) 
Fit Orig. 


Ml' 2 (10- 3 /m.+ ) 
Fit Orig. 


0.20 SAID 
SL 
DMT 
MAID 


-1.41 ±0.67 -1.47 
-2.24 ±0.69 -3.11 
-1.75 ±0.67 -2.82 
-2.36 ±0.69 -2.06 


-4.68 ±0.28 -5.85 
-5.11 ±0.27 -4.64 
-5.04 ±0.27 -6.65 
-5.50 ±0.29 -6.50 


38.89 ± 0.44 ± 0.62 39.85 
39.76 ± 0.43 ± 0.62 40.48 
39.84 ± 0.43 ± 0.62 39.65 
39.43 ± 0.43 ± 0.62 39.98 


Avg. 


-1.96 ±0.68 ±0.01 ±0.41 


-5.09 ±0.28 ±0.02 ±0.30 


39.57 ±0.43 ±0.62 ±0.40 


GH 
PV 


-1.68 ±0.47 
-3.05 ± 1.20 


-6.75 ± 1.85 
-9.19 ±3.00 


38.22 ± 5.10 



TABLE III: Values of EMR, CMR, and M 1+ at W = 1232 MeV and Q 2 = 0.060 (GeV/c) 2 (top) and Q 2 = 0.200 (GeV/c) 2 (bot- 
tom) for the EFT predictions and fitted models. The uncertainties are in the order of statistical then systematic. The systematic 
uncertainties for the individual models' EMR and CMR are suppressed because they are small. For the average, the third num- 
ber is the model uncertainty defined as the RMS deviation of the results from the four different models. The models are the 
three resonant parameter fitted SAID [lj, MAID [37[, Sato-Lee(SL) and DMT [H, [HI] models at W = 1232 MeV (1227.3 
MeV for SAID [13]) and Q 2 = 0.060 (GeV/c) 2 . The EFT predictions are Gail and Hemmert (GH) [H and Pascalutsa and 
Vanderhaeghen (PV) [13, [H|| and are presented without fitting to data. 
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Parallel cross section 



In Figures [TTI and [T2l the parallel cross section W scans 
at Q 2 = 0.060 and 0.200 (GcV/c) 2 , respectively, are plot- 
ted along with corresponding model predictions. In Panel 

(a) of Fig. [11] the unmodified models are shown. In Panel 

(b) of Fig. [TTJ the results of the three resonant param- 
eter fit to the previously shown data were used. It is 
important to note that in Fig. [Til only the helicity inde- 
pendent, low Q 2 results from Sec. IV Al have been fit and 
yet the agreement with the W scan data is improved sig- 
nificantly. However, there is still disagreement with the 
data near the tails which indicate issues with the model 
backgrounds. 

The model curves in Fig. [T2] were made in the same 
way as those for Fig. [TTJ Again, even though the W 
dependent data were not included in the fit, the models 
converged noticeably. There are even larger deviations 
at high W indicating an additional Q 2 dependence to 
the model background terms which is not accounted for 
properly. It is hoped that both these sets of W dependent 
data will help to constrain the models once the models 
have been improved. 

One property that these data can help determine is the 
shape of the parallel cross section versus Q 2 in the range 
from Q 2 = 0.060 to 0.200 (GeV/c) 2 . The four models 
(MAID, Sato-Lee, SAID, and DMT) do not have a large 
variation of the shape of the parallel cross section with 
Q 2 but they do differ from one another in peak center 
value and width. This was found by plotting the model 
predictions versus W for the Q 2 range of the data after 
normalizing the cross sections to the value at the peak. 
The same procedure was carried out on the Mainz and 
Bates data and plotted with the peak normalized DMT 
model in Fig. 1131 To aid in the comparison, fits to the 
data were performed using a Breit-Wigner form plus a 
quadratic background. The fits determined the peak to 
be W = 1206 ± 1 MeV with widths of 83 to 108 MeV 
with 10% uncertainties. The shape of the data clearly 
does not change dramatically. There does appear to be 
some deviation from the predicted shape toward high W 
in this model. 

In addition to the data taken at Q 2 = 0.060, 0.127, and 
0.200 (GeV/c) 2 , one data point was taken at Q 2 = 0.300 
(GeV/c) 2 . This was a parallel cross section measurement 
near the peak of the cross section, W — 1205 MeV. As 
can be seen in Fig. [TJJ all of the points taken at Mainz 
are consistent with the unfit Sato-Lee model but show 
the same variation with Q 2 as all the models. Previous 
data from Bates [l5[ tend to have larger values in general 
as in Fig. [7] but are within the combined uncertainties. 
The variation with Q 2 is significant because the shape 
is consistent with a large pion cloud contribution. Note 
that none of the models in Fig. [14] were fit to the data 
and the spread in their predicted values is similar to the 
spread seen in the Panels (a-d) of Figs. [Til and [T2l which 
also show unfit models. 
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FIG. 11: (Color online) Parallel cross section for the 
p(e, e'p)-K° reaction at Q 2 = 0.060 (GeV/c) 2 before (panel 
a) and after (panel b) three resonant parameter fit. Model 
curves are the same as in Fig. [8] The smaller error bars are 
the statistical uncertainty and the larger error bars include 
the systematic uncertainty added in quadrature. 
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FIG. 13: (Color online) Parallel cross section data (uncertain- 
ties are statistical and systematic added in quadrature) for the 
p(e, e'p)iv reaction from Mainz (this paper) and Bates [HI, [13] 
and model predictions from [2j| scaled so that the W = 1205 
MeV peak is at 1.0. Note that the shape of the parallel cross 
does not vary by a large amount. The data points are con- 
nected by fits (thick lines) which were a Breit-Wigner form 
plus a quadratic background. 
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FIG. 12: (Color online) Parallel cross section for the 
p(e,e'p)ir° reaction at Q 2 = 0.200 (GeV/c) 2 before (panel 
a) and after (panel b) fit. Model curves are the same as in 
Fig. [8] The uncertainty is the statistical and systematic un- 
certainties added in quadrature. 
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FIG. 14: (Color online) Results for cr versus Q 2 for W = 1205 
MeV, 6p q = 0°. The solid circle data were taken at Mainz and 
include the statistical and systematic uncertainties added in 
quadrature. The Q 2 = 0.127 (GeV/c)f_point was taken at 
W = 1212 MeV and models g3, HI, H3, HJ were used to 
extrapolate to W = 1205 MeV. As a result, the uncertainty 
is slightly larger for that point than the others. The plotted 
models have not been fit to the data. 
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C. Background Sensitive Data Below Resonance 

Background sensitive data were taken at Q 2 = 0.060 
(GeV/c) 2 and at low W, where the M 1+ amplitude is 
less dominant, to test the background amplitudes in the 
reaction models. Comparing over a wide range of W is a 
rigorous test of the background multipoles and the shape 
of Mi+. Also, the background multipoles are more im- 
portant at low W where they are relatively larger than 
the resonant multipoles which are then off-resonance. In 
addition, the -A/i+ term is not purely imaginary in that 
region and interferences from real background amplitudes 
will not be suppressed as much. Although the extrac- 
tion of specific background multipoles was not planned, 
model predictions can be compared to the data to see 
whether they agree or not. In addition, some fitting in- 
cluding background terms can provide an indication as 
to which amplitudes may be significant. Similar stud- 
ies were performed in Ref. [13] on the Bates Q 2 — 0.127 
(GeV/c) 2 data. 

Figure [15] shows the new low W data at Q 2 = 0.060 
(GeV/c) 2 compared with models before and after fitting 
the three resonant multipoles at resonance. In the un- 
fit panel (a) of Fig. [15] only the DMT model is close 
to the data but does not reproduce both points. The 
other models differ in size and shape due to their differ- 
ent background amplitudes. The same low W data points 
are plotted in the lower panel of Fig. [15] along with the 
model fits at resonance found in Sec. IV Al These data 
points were not included in those fits. The agreement 
with the models has not improved significantly. This is 
corroborating the known disagreement among the models 
for the various background amplitudes. 

Taken together, both panels of Fig. [IS] are indicating 
that, as expected, more than a three parameter resonant 
fit is required away from the resonance region. As an 
example, one of the best four parameter fits (using the 
three resonant parameters along with which in- 

clude these background data is shown in the lower panel 
of Fig. [151 While that fit is satisfactory for the two low 
W data points, it is still in disagreement with the back- 
ground sensitive ctlt' data near resonance and for some 
of the parallel cross section W scan results especially at 
the low and high W tails. Because of the sensitivity of 
these few points, the overall x 2 was not improved very 
much. However, as mentioned above, since the data set 
is limited, it is not a surprise that no single background 
multipolc allows a good fit. The effect of background am- 
plitudes from the models can be compared to data but 
the amplitudes themselves cannot be determined. 

Figure fTS] indicates the need for more precise model 
calculations and possible estimates of uncertainties. In 
addition, dedicated low W experiments could help con- 
strain the models and lead to more refined predictions. 

After comparing the data and models over a range of 
observables, the DMT model has the best overall agree- 
ment with all of the low Q 2 data. The fitted DMT result 
in Fig. [TSJis the closest to the data of all the models. The 
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FIG. 15: (Color online) Background sensitive data from 
Mainz at W = 1155 MeV, Q 2 = 0.060 (GeV/c) 2 before 
(panel a) and after (panel b) fit to the data near resonance. 
The smaller error bars are the statistical uncertainty and the 
larger error bars include the systematic uncertainty added in 
quadrature. The curves are MAID 2003 (dotted) [13 DMT 
(dot-dash) [Ilfjsato-Lee (dot-dot-dash) [27| and SAID (dot- 
dot-dot-dash) [5(J. Panel (b) also includes an example of one 
of the four parameter fits (long dash) including these back- 
ground data. That fit used the three resonant parameters 
along with Mi-) using the MAID 2003 model. 



fitted DMT results for o^T' hi Fig. [5] are fairly close to 
the data and no worse than MAID and SAID. Finally, the 
fitted DMT results for the W scan in Fig. [11] look very 
good overall and only disagree at a few points. While 
no model agrees perfectly with all the data, the DMT 
model after the three resonant parameter fit does appear 
to describe the new data the best. 

Studies using various fitting parameters indicate a path 
to follow for improving the agreement between data and 
theory. The problem encountered with the fitting method 
used in this work is that the parameters apparently do 
not give the models enough freedom to fit the background 
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amplitudes. As is clear in the plot in Fig. [ITJ the models 
simply have the wrong shape. The next step to con- 
strain the background is to fit the less well determined 
coupling constants and other internal model terms which 
affect many multipoles at once. This would hopefully 
add enough freedom to allow the models to fit the data. 
Despite these quantitative problems, the background am- 
plitudes are sufficiently small near resonance so that the 
uncertainties in them do not contribute more than the 
experimental uncertainties in determining the resonant 
amplitudes. 



VI. Q 2 VARIATION OF RESONANT 
MULTIPOLES 

One of the main goals of this experiment is to deter- 
mine the Q 2 variation of the resonant multipoles. The 
new Jefferson Lab results [(| for the EMR and CMR 
at Q 2 = 0.200 (GeV/c) 2 agree with our results very 
well. In addition, it was already shown that the present 
Mainz data agree with the previous Bates data [l5| at 
Q 2 = 0.127 (GeV/c) 2 . All of this shows that there is 
reasonable consistency of the results from the different 
laboratories. 

Figure [16] shows the evolution of the multipoles at 
low Q 2 along with all the other published points. Two 
representative constituent quark models, the newer hy- 
percentral quark model (HQM) [6(3|, and an older non- 
relativistic calculation of Capstick [23[, have been in- 
cluded (the relativistic calculations are in even worse 
agreement with experiment). These curves are represen- 
tative of quark models which typically undcr-prcdict the 

dominant multipole by ~ 30% and underestimate 

the EMR and CMR by an order of magnitude, even pre- 
dicting the wrong sign. One solution to this problem 
has been to add pionic degrees of freedom to quark mod- 
els [MUS El- All of these models treat the A as a bound 
state and therefore do not have the 7rN continuum (i.e. 
no background amplitudes ) so that cross sections are not 
calculated. The Sato-Lee [27| and DMT [H] dynamical 
reaction models with pion cloud effects bridge this gap 
and arc in qualitative agreement with the Q 2 evolution of 
the data. These models calculate the virtual pion cloud 
contribution dynamically but have an empirical param- 
eterization of the inner (quark) core contribution which 
gives them some flexibility in these observables. By con- 
trast the empirical MAID [37| and SAID [5(| represent 
fits to other data with a smooth Q 2 dependence. 

Both the dynamical [2?], [H[ and the phenomenolog- 
ical [37], H3] models are in qualitative agreement with 
the experimental results. Nevertheless, all models exhibit 
some small deficiencies either on top or at the wings of the 
resonance indicating that detailed improvements could 
and should be implemented to the models description 
of resonant or background amplitudes towards account- 
ing for these deficiencies. As a general remark one can 
note the much better behavior of the dynamical models 



(DMT and Sato-Lee) compared to the phenomenological 
ones (MAID and SAID) as far as the description of the W 
evolution of the cross section is concerned (see Fig. ITTj) 
while for Sato-Lee the description of the fifth response 
is also excellent thus indicating that the model provides 
the most consistent description of the background am- 
plitudes. One must also point out though the consistent 
description that SAID provides for the unpolarized cross 
sections on top of the resonance measurements for all Q 2 
points. 

The plotted lattice QCD results with a linear pion mass 
extrapolation (58j are in general agreement with the data 
for the EMR but disagree for the CMR by a wide mar- 
gin. This margin is bridged, though, when using a chiral 
extrapolation to the physical pion mass instead of the 
linear one. The EFT analysis of Pascalutsa and Vander- 
hacghen (PV) [IcJ indicates that a linear extrapolation 
is close to the data for the EMR but not for the CMR for 
which these extrapolated lattice results are considerably 
reduced. The second plotted lattice QCD results were 
performed with an improved method an a smaller pion 
mass and are reported without any extrapolation [59( . It 
is significant that these newer results have the same sign 
as the data at low Q 2 . The general qualitative agreement 
of the lattice QCD calculation provides a direct link with 
the experimental evidence for deformation to QCD. 

The results of the two effective field theory calcula- 
tions [29|, [3(| are also presented in Fig. [T6j These con- 
tain empirical low energy constants. For Gail and Hcm- 

3/2 

mcrt this includes fits to the dominant M 1+ multipole 
for Q 2 < 0.2 (GeV/c) 2 and for the EMR at the pho- 
ton point (Q 2 =0). In order to achieve the good over- 
all agreement they had to employ one higher order term 
with another empirical constant. The EFT calculation 
of Pascalutsa and Vanderhaeghen [3(j provided a valu- 
able estimate of the uncertainties caused by excluding 
the next higher order terms from the calculation. While 
this is a very helpful start, the uncertainties are signif- 
icantly larger than the experimental uncertainties and 
will have to be reduced through a proper treatment of 
the excluded higher order terms. However, these effec- 
tive field theoretical (chiral) calculations that are solidly 
based on QCD, successfully account for the magnitude 
of the effects giving further credence to the dominance of 
the meson cloud effect. 

One way to see the major role played by the pion cloud 
contribution to the resonant multipoles is that for this 
case the expected scale for the Q 2 evolution is m 2 = 0.02 
GeV 2 . In these units the range of the present experi- 
ment for Q 2 from 0.060 to 0.200 (GeV/c) 2 is 3 to 10 
units. Therefore it is not surprising that one should sec 
relatively large changes in the predicted Q 2 evolution of 
the resonant multipoles as is shown in Fig. [T6j It is also 
clear that there is significant model dependence in these 
predictions. 
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FIG. 16: (Color online) The low Q 2 dependence of the Mi+, EMR, and OMR at W = 1232 MeV for the 7*p -» A reaction. 
The T symbols are our data points and include the experimental and model uncertainties (see Table [TTT)l added in quadrature. 
The other data are the photon point data O [3 and (g> [Ic|, CLAS □ [1 EH (CLAS data for 0.16 < Q 2 < 0.36 (GeV/c) 2 are 
from a unitary isobar model fit and have statistical uncertainties only), Bates A [13], Eisner ® [H}> an d Pospischil EE 
All uncertainties are statistical and systematic added in quadrature unless otherwise noted. The lattice QCD calculations with 
linear pion mass extrapolations are shown as x [k| and the new calculations with small pion mass but without extrapolation 
are shown as • [kJ. Also shown are the chiral perturbation calculations of Pascalutsa and Vanderhaeghen (PV) (see EFT in 
Fig. [8} and Gail and Hemmert (GH) (block solid line) [2iJ. The other curves represent the same models as in Fig. [8] The 
HQM (long-dashed line) [60l ] and Capstick (short-dashed line) [13] quark models have been included. 



VII. CONCLUSIONS 

The data presented here provide a precise determina- 
tion of the resonant amplitudes in the j*p — > A reaction 
in the range of Q 2 = 0.06 to Q 2 = 0.20 (GeV/c) 2 (3 to 
10 m\). The experiment at the Mainz Microtron was 
carefully designed to reach the lowest possible Q 2 to test 
effective field theory calculations and to probe the regime 
where pionic effects are predicted to be a maximum and 
to vary significantly [U HI] . The absolute cross section 
accuracy at the 3% level was verified with several cross- 
checks. The measurement of the <tq = <tt + £<tt, <Jlt, &tt 
and <jlt> partial cross sections, at center of mass ener- 
gies both on and off resonance, allows for sensitive tests of 



effective field theory pjj , l30l | and reaction model calcula- 
tions [13, [H, H3, Hoj ■ These partial cross sections are also 
important for extracting the resonant multipoles from the 
data; these are used to test lattice calculations [H, |59j] 
and quark models [6(| [HI • At the present time the ex- 
periments are more accurate than both theory and model 
calculations. 

The chiral effective field theory predictions [1^, [3C| 
agree with our cross section data within the relatively 
large estimated theoretical uncertainties due to the ne- 
glect of higher order terms. It is clear that a quanti- 
tative comparison of these calculations and experiment 
must wait until the next order calculations are per- 
formed. The phcnomenologically adjusted models Sato- 
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Lee, DMT, SAID and MAID [13, H, 13, Hq| are in 
good agreement with experiment when the resonant am- 
plitudes are adjusted to the data. This allows an accu- 
rate extraction of the Ml, E2 and C2 resonant multipolcs 

3/2 3/2 3/2 

(M 1+ , E 1+ , S 1+ ) with an estimated model uncertainty 
which is approximately the same as the experimental un- 
certainty. This has been achieved due to the precision of 
the experimental data and also because of the dominance 
of the magnetic dipole amplitude A/ 1+ ; this dominance 
means that differences in the background amplitudes are 
not significant near resonance and that the model uncer- 
tainties in the determination of the resonant multipolcs 
are comparable with with the experimental uncertain- 
ties. The differences in the background amplitudes have 
been demonstrated in our low W data and in <7lt' for 
which the background multipolcs play a more significant 
role. This emphasizes the need for model builders to im- 
prove their calculations and also to present their uncer- 
tainties, as has been done in the EFT calculations. Wc 
have performed our own error estimate by comparing the 
extracted resonance multipoles using different models. 

Comparisons of the measured resonant multipoles as a 
function of Q 2 show reasonable agreement between ex- 
periments at different laboratories. The non-zero values 

3/2 3/2 

of the quadrupole amplitudes (E 1+ , S 1+ ) demonstrate 
the existence of non-spherical amplitudes in the nucleon 
and A conjectured many years ago on the basis of the 
non-spherical interaction between quarks plj . This fea- 
ture is also present in the lattice calculations [H, \E^ , 
thus linking the experimental evidence for deformation 
directly to QCD. Unfortunately, the uncertainties in the 
present calculations are large, which precludes a quanti- 
tative comparison with experiment. We anticipate fur- 
ther advances with calculations at lower quark masses 
combined with improved chiral calculations which are 
also just in the beginning [30| . These results show quali- 
tative agreement with the two chiral effective field theory 
results [H, 30] . The uncertainties in these latter two cal- 
culations indicate that higher order terms must be eval- 
uated before a quantitative comparison can be made. 

Comparison with representative quark models [60L l6l| 
shows that they are not close to the data indicating 
a deficiency of the underlying physics description while 
demonstrating that the color hyperfine interaction is in- 
adequate to explain the effect, at least at large dis- 
tances. Our present understanding is that the long range 
(low Q 2 ) region is dominated by the spontaneous break- 
ing of chiral symmetry in QCD which results in non- 
spherical pion emission and absorption from the nucleon 
and A [3,1,111. 

Even though experiments are ahead of theory at the 
present time, future experiments can add to the current 
understanding by measuring in the mr + and 7 channels 
and by utilizing polarized targets and polarimeters. New 
data for the 7 channel from Mainz are under analysis 
and should be published soon [(52|. Additional low W 
data would also give a better handle on the background 
amplitudes. 
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APPENDIX: DATA TABLES 
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TABLE IV: Q 2 = 0.060 (GeV/c) 2 cross sections. The first 
uncertainty is the statistical uncertainty and the second is 
the systematic uncertainty. The helicity dependent cross sec- 
tions are shown in Table FVl The uncertainties are statistically 
dominated for those results. The third set of results are from 
the background amplitude test. The lower set of results is the 
W parallel cross section scan. 
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TABLE V: Summary of the extracted values for <7o, ffrr, 
(tlt, and a LT i. The uncertainties in the cross section are the 
statistical and systematic uncertainty respectively. See text 
for details of the uncertainty estimation procedure. 
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TABLE VIII: Q 
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0.200 (GeV/c 


I 2 cross section results. The 



uncertainties correspond to the statistical and the systematic 
uncertainties respectively. The lower set of results is the W 
parallel cross section scan. 



TABLE VI: Q 2 = 0.127 (GeV/c) 2 results. The first uncer- 
tainty is statistical and the second is the systematic. The 
helicity dependent results are in Table [VTT1 
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TABLE VII: Summary of the extracted values for <ro, olt, 
and a LT i at Q 2 = 0.127 (GeV/c) 2 . The uncertainties in the 
cross section are the statistical and systematic uncertainty 
respectively. 
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TABLE IX: Extracted values for oo, ott, olt, and o L t' a t 
Q 2 = 0.200 (GeV/c) 2 . The uncertainties correspond to the 
statistical and the systematic uncertainties respectively. 
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TABLE X: Q 2 = 0.300 (GeV/c) 2 results. The first uncer- 
tainty is statistical and the second is systematic. 
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